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ABSTRACT 

We have recently carried out the first wide-field hydrogen Paschen-a line imaging 
survey of the Galactic Center (GC), using the NICMOS instrument aboard the Hubble 
Space Telescope. The survey maps out a region of 2253 pc 2 (416 arcmin 2 ) around 
the central supermassive black hole (Sgr A*) in the 1.87 and 1.90 urn narrow bands 
with a spatial resolution of ~ 0.01 pc (0'.'2 FWHM) at a distance of 8 kpc. Here 
we present an overview of the observations, data reduction, preliminary results, and 
potential scientific implications, as well as a description of the rationale and design 
of the survey. We have produced mosaic maps of the Paschen-a line and continuum 
emission, giving an unprecedentedly high resolution and high sensitivity panoramic 
view of stars and photo-ionized gas in the nuclear environment of the Galaxy. We detect 
a significant number of previously undetected stars with Paschen-a in emission. They 
are most likely massive stars with strong winds, as confirmed by our initial follow- 
up spectroscopic observations. About half of the newly detected massive stars are 
found outside the known clusters (Arches, Quintuplet, and Central). Many previously 
known diffuse thermal features are now resolved into arrays of intriguingly fine linear 
filaments indicating a profound role of magnetic fields in sculpting the gas. The bright 
spiral-like Paschen-a emission around Sgr A* is seen to be well confined within the 
known dusty torus. In the directions roughly perpendicular to it, we further detect 
faint, diffuse Paschen-a emission features, which, like earlier radio images, suggest an 
outflow from the structure. In addition, we detect various compact Paschen-a nebulae, 
probably tracing the accretion and/or ejection of stars at various evolutionary stages. 
Multi-wavelength comparisons together with follow-up observations are helping us to 
address such questions as where and how massive stars form, how stellar clusters are 
disrupted, how massive stars shape and heat the surrounding medium, how various 
phases of this medium are interspersed, and how the supermassive black hole interacts 
with its environment. 

Key words: Galaxy: center; Interstellar Medium (ISM); infrared: stars; circumstellar 
matter; stars: formation 



1 INTRODUCTION 

The Galactic center (GC) is a unique laboratory for a de- 
tailed study of star formation and its impact on the nuclear 
environment of galaxies. Because of its proximity, the GC 
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can be imaged at resolutions more than 100 times better 
than the nearest galaxies like our own (e.g., the Andromeda 
galaxy). This unmatched high-resolution capability is par- 
ticularly important to the study of the formation and de- 
struction processes of stellar clusters as well as the overall 
population of massive stars in such an environment. Indeed, 
the GC is known to host three very massive clusters (the 



2 



Arches and Quintuplet clusters, as well as the central par- 
sec cluster surrounding Sgr A*; ages ~ 2 — 7x 10 6 yrs, masses 
~ 3xl0 4 M e each; e.g. ,|Cotera et al.|l996[|Serabyn, Shupe fc 
Figer|1998l|Figer et al.|1999l|2002||Genzel et al.|2003h, which 



are responsible for about 5% of the total Lyman continuum 
flux for the Galaxy, in less than 0.01% of its volume. Mas- 



sive stars are known to exist outside the clusters ( Cotera et 
|al.||2006j ), however, locating these important stars has been 



sporadic and incomplete (e.g. Muno et al. 2006a Mauer 



han, Muno & Morris 2007 Mauerhan et al. 20091. These 



isolated massive stars may have formed individually or in 
small groups, although the extreme conditions in the GC 
make this mode of star formation unusually difficult. They 
may also have been spun out of the clusters due to their in- 
ternal dynamics. This process is expected to be accelerated 



by the strong external tidal force in the GC (Kim, Morris & 
Lee|1999p , 



The star formation process depends critically on the in- 
terplay of massive stars with the interstellar medium (ISM). 
Massive stars release large amounts of ultraviolet radiation. 
They also produce strong stellar winds throughout their 
short lives and die in supernova explosions. Such energetic 
stellar output dramatically shapes the morphology of the 
ISM, as is revealed in existing radio, mid-infrared, and X- 
ray maps (Fig. [T| e.g 



Yusef-Zadeh, Morris & Chance 



T984l|Yusef-Zadeh et al.|2002||Wang, Gotthelf fc Lang|2002 



Wang, Dong & Lang||2006 Muno et al. 112009 \. Perhaps the 



best-known example of this is the Radio Arc region which 
surrounds the Arches and Quintuplet clusters: a collection of 
arc-like thermal filaments, synchrotron-bright non-thermal 
linear filaments, and diffuse, hot gas. Some of the structures 
within this region, the Sickle HII region ionization front and 
photoevaporated pillars in particular, suggest that gas is be- 
ing collected and compressed ( Cotera et al.|2006 Cotera et 
al. in preparation), and is morphologically resembles known 
regions outside the GC where such processes are believed 
to have led to the formation of a new generation of stars 
(e.g. M16;|Smith, Stassun fc Bally| ([20051)). The ener gy in- 



jection from massive stars and possibly from the central su- 
permassive black hole (SMBH), that is coincident with the 
radio source Sgr A*, may also strongly affects the thermal 
and/or dynamical properties of the ISM. Molecular gas in 
the GC, for example, has unusually high turbulent velocities 
and temperatures compared to clouds in other parts of the 



Galaxy (e.g., Morris & Serabyn 19961. This may bias the 
stellar initial mass function toward heavier stars and favor 
the production of massive star clusters ( Morris [1993 Wang, 
Dong fc Lang|2006| and references therein). 

To constrain the dynamical process of the known clus- 
ters, the overall population of massive stars, and their forma- 
tion modes, we sought to obtain a uniform survey of massive 
stars and their interplay with the ISM across the central 
90 parsecs. Ionized warm gas provides the ideal tracer of 
massive stars and their interplay with their surroundings. 
Although previous radio continuum observations are free of 
extinction and trace both thermal and non-thermal emis- 
sion, the scale-dependence of radio interferometers makes 
it difficult to accurately measure flux densities on arcsec- 
ond scales. In addition, fine features (e.g., ionization fronts) 
are often found in the midst of diffuse non-thermal emis- 
sion (e.g., the Radio Arc), making decomposition of thermal 



and non-thermal components challenging. Radio recombina- 
tion line observations typically have a substantially lower 
signal-to-noise and spatial resolution than radio continuum 
data. In the near-IR, the 2.16^m Br7 line is accessible from 
the ground (no suitable filter is available for existing and 
planned space-based platforms), but has also not proved ef- 
fective for large scale mapping of warm ionized gas. To our 
knowledge, no large-scale Br7 line survey of the GC has ever 
been published. The difficulty lies in the lack of the com- 
bination of high spatial resolution and stable point spread 
function (PSF) in a ground-based survey. Observations with 
adaptive optics (AO) can have superb resolution, but only 
over small fields of view (typically 15" - 40"). The AO PSF 
also changes strongly, both spatially and with time, and has 
extended PSF wings, making the PSF subtraction extremely 
difficult. In addition, photometric conditions rarely, if ever, 
persist for the time required for such a survey, hampering 
the study of large-scale features. 

We have carried out the first large-scale, high- 
resolution, near-IR survey of the GC, using HST NICMOS. 
The primary objective of the survey is to map out the hydro- 
gen Paschen-a (Pa for short) line (wavelength 1.876 \im) in 
a field of ~ 39' x 15' around Sgr A* (corresponding to 90 pc x 
35 pc at the GC distance of 8 kpc; |Ghez et al.|2008| Fig.[l]). 
This field is known to be rich in clustered massive star for- 
mation (i.e., including the three known clusters). NICMOS 
observations of the Pa emission were taken previously only 
for a few discrete regions, centered on the three known mas- 
sive clusters (see § 4.3 for further discussion). Our survey 



field also includes much unexplored territory. In particular, 
the region between Sgr A and Sgr C along the Galactic plane 
(to the right in Fig. [TJ shows a number of compact Spitzer 
IRAC sources (several with radio counterparts) likely related 
to some of the earliest stages of massive star activities. This 
field selection also facilitates an unbiased study of the appar- 
ent asymmetry between the positive and negative Galactic 
longitude sides relative to Sgr A*. 



Our HST NICMOS observations overcome all of the po- 
tential problems of a ground-based survey. The wide NIC3 
camera on this instrument provides a resolution of ~ 0"2 
with a stable PSF. This capability is important for detecting 
and cleanly removing point-like sources as well as resolving 
fine structures of extended features. Both the stable PSF 
and the absence of atmosphere for our HST survey means 
that the entire survey is photometrically accurate and con- 
sistent. The extra magnitude of extinction along the sight- 
line to the GC at 1.87u.m (as compared to 2.16 um) is more 
than compensated for by an intrinsic Pa line intensity 12 
times greater than the Br7 line ( Hummer fc Storey||1987 |. 
In addition, the background at Pa with NICMOS is a factor 
of ~ 800 lower than the sky background at Br7 as observed 
from the ground. This low background of HST NICMOS also 
allows us to detect faint point-like sources and extended fea- 
tures, particularly over our large mosaicked field. 



2 OBSERVATIONS AND DATA REDUCTION 

Our survey was completed in 144 HST orbits between Febru- 
ary 22 and June 5, 2008. All of these orbits had an identi- 
cal design, each consisting of 2 x 2 NIC3 pointing positions 
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Figure 1. A multi-wavelength montage of the GC: the VLA 20 cm continuum (red; [Yuscf-Zadch, Morris & Chance 1984), Spitzer 8 
jim ( green; |Stolovy et"aT1|2006| |Arendt et aL]|2008t , and Chandra ACIS-I 1-9 keV (blue; |Wang, Gotthelf fc Lang||2002| |Muno et al.| 
|2009) . The three known massive young clusters are shown, as well as the outer border of the Pa survey. 



for both F187N (on-line) and F190N (off-line) filters. These 
filters have a 1% bandpass centered at the wavelengths 
1.875um and 1.900 \im. NIC3 is a 256x 256 HgCdTe array 
detector with a pixel size of ~ 0'.'2031 x 0'.'2026 in the x and 
y directions. Adjacent positions and orbits were designed to 
have 3" overlaps in the vertical direction of the NIC3 array 
coordinates, partly to allow for the bottom 15 pixels that 
do not provide useful data due to afocal vignetting of that 
portion of the array. This is in addition to the overlaps re- 
sulting from the 4-point dithers at each of the four positions 
of an orbit. Each dither exposure was taken in the MULTI- 
ACCUM readout mode (SAMP-SEQ=16; NSAMP=9) for 
photometry over a large flux range and for effective iden- 
tification and removal of cosmic-rays and saturated pixels. 
The 4-point dithers were in an inclined square wave pat- 
tern (i.e., the shifts were in the directions of 67°. 5 relative 
to the x-axis of the detector coordinates); each dither step 
produced a shift of 2'.'33 and 5'.'63 pixels in the x and y di- 
rections. This dither pattern enabled both a sub-pixel sam- 
pling of the undersampled NIC3 PSF and a consistent over- 
lap between orbits. Each orbit effectively covered an area of 
roughly 102" x 98" (2.8 arcmin 2 ). With this strategy and a 
single telescope orientation, the survey mapped out a field 
of 416 arcmin 2 (as outlined in Fig. [I]), and with a uniform 
192 sec exposure per filter. The entire science data set in- 
cludes 144 (orbits) x 4 (pointing positions) x 4 (dithers) x 
2 (filters) = 4604 images. 

In addition, between target visibility periods (during 
Earth occultations) after each of the 144 orbits, 16 dark 
frames were obtained with the detector clocked in a manner 
identical to the science exposures. These inter-orbit raw dark 
exposures were assembled to form a "superdark" frame used 
in our image calibration. 

The data reduction and analysis, based on a combi- 
nation of several newly developed tools and the standard 
routines implemented in the software IRAF/STSDAS, will 



standard DC offset correction, we use a global fitting method 
to determine the relative instrumental background offsets 
among all NIC3 array quadrants, pointing positions, and or- 
bits, based on multiple overlapping fields. A similar method 
is also used to optimize the correction for the relative as- 
trometry offsets among orbits. The absolute astrometry is 
corrected in respect to radio counterparts of the near-IR 
sources in the field, reaching an accuracy better than 0.1". 
The source detection is based primarily on the Starfinder 
software ( jDiolaiti et al.|2000[ ). 

The procedure and product of the data reduction are 
illustrated in Fig. [2] for global views and Fig. [3] for full- 
resolution close-ups. To estimate the continuum contribu- 
tion in the F187N band (Figs. [2^, and [3^), we adaptively 
construct a mean F187N to F190N intensity ratio map, with 
the value of each pixel given as the median ratio of the near- 
est 100 sources, most of which are assumed to be located 
at the GC. However, this assumption does not hold for re- 
gions with exceptionally high extinction, as judged from a 
low surface number density of detected sources (i.e., <0.3 
sources arcsec -2 ). For such regions, we calculate the ratio 
values from the extinction map obtained from Spitzer obser- 
vations ( |Schultheis et al.|2009||Dong et al |2009| for details). 
The continuum-subtracted map (Figs. [2jj and [SJa) gives a 
crude look of the potential Pa-emitting sources as well as the 
diffuse emission. However, one cannot judge whether or not 
a peak (e.g., in Fig.JSJj or]!]}) is a Pa emission source based 
only on its apparent excess above the background; such a 
peak can easily be caused by the statistical flux fluctuation 
of a subtracted source, especially a bright one. In fact, all 
peaks in[3j3, except for the brightest one, are below our cur- 
rent Pa source detection threshold. We identify Pa source 
candidates based on their integrated continuum-subtracted 
fluxes at a statistical confidence > 5a, where a accounts for 



be detailed in Dong et al. (20091. Briefly, in addition to the 



all known measurement errors (Dong et al. 20091. To pro- 
duce a diffuse Pa map (Figs. |2]: and |3b), we use the ac- 
tual F187N to F190N flux ratio of a detected source for 
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Figure 2. NICMOS maps from our survey: (a) calibrated F187N mosaic, (b) continuum-subtracted F187N map, and (c) continuum- 
subtracted and point source removed F187N (diffuse Pa) map. These maps are presented with a bin size of 0''4. The intensity is 
logarithmically scaled from of 0.01 mjy arcsec -2 to 10 mjy arcsec -2 in (a) and to 3 mjy arcsec -2 in (b) and (c). 
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its affected pixels (defined to have its PSF-predicted inten- 
sity contribution > 5cr above the local background in the 
F190N map). We further remove any uniformly-distributed 
foreground and residual instrument background in a map 
by shifting its intensity to zero, averaged over a few darkest 
regions (see §[3). 



3 PRELIMINARY RESULTS 

Fig. [2] shows a panoramic view of our survey data. Clearly, 
the distribution of both the stellar light and the diffuse 
emission is highly inhomogeneous, partly due to the vari- 
ations in foreground extinction, which varies strongly from 
one region to another and on different spatial scales. The 
most outstanding extinction features in Fig. [2] are the dis- 
tinct dark filaments (tendrils) in the field just above Sgr A, 
which are also noticeable in published K-band images (e.g., 
Philipp et al. 1999[ ). The darkness of these filaments indi- 
cates that they are foreground thick dusty clouds silhouetted 
against the IR-luminous nuclear region of the Galaxy The 
filaments appear to be structured down to arcsecond scales 



Figure 3. A full-resolution (O'.'l pixel) close-up of a newly dis- 
covered Pa-emitting complex (G359. 866+0. 002), about 5' north- 
west of Sgr A*: (a) F187N image, (b) continuum-subtracted 
F187N image, and (c) point source removed (diffuse Pa) image. 
The intensity is logarithmically scaled over the range of 0.1-30 
mjy arcsec -2 in (a) and 0.1-1 mjy arcsec -2 in (b) and (c). 



[or ~ 10 1 (d/lkpc) cm, where d is the distance to such a fil- 
ament; see the full resolution close-ups in Fig. [3]. With our 
data, we can distinguish between extinction and lack of ion- 
ized emission by comparing the Pa and F190N images - if 
a region is dark in both images, we can say with confidence 
that the dark features indicate regions with extinctions Av 
> 40 - 50. 



3.1 Extended Pa Emission 

Fig. [2ji shows extended Pa-emitting complexes across 
much of the field. Strong and distinctive emission fea- 
tures are seen on angular scales from sub-arcsecond to 
about 10', with surface brightness ranging from as high 
as ~ 10 mjy (millijansky) arcsec -2 down to a 3<r level of 
~ 0.05 mjy arcsec -2 for the image with the 0'.'4 x 0'.'4 
pixel. The rms noise of a map depends on its pixel size and 
is ~ 0.06 mjy arcsec -2 for the full resolution image with 
a pixel size of O'.'l x O'.'l, accounting for both statistical 
and systematic fluctuations (photon counting, artifacts from 
stellar continuum subtraction, etc.). This noise is estimated 
from selected regions in the Pa map (e.g., at the upper 
right corner of the map; see Dong et al.||2009 for details) 
that show little Pa emission but are bright in the F190N 
continuum image (hence the rms contribution from clumpy 
structures in the Pa emission or extinction should be min- 
imal). The overall distribution of the Pa emission is highly 
lopsided toward the positive Galactic longitude side of the 
GC, similar to what is seen in radio and mid-IR (Fig. [T]). 

The diffuse Pa emission is greatly enhanced around 
the three known clusters responsible for the ionization of 
the inner rims of surrounding gas. The emission around the 
GC cluster is largely confined within a small spiral-like nu- 
clear structure, the Pa emission of which was first detected 
by ( |Scoville et al.|2003[ ). Our new survey now enables us to 
trace the faint emission features to large distances, which 
extend vertically below and above this structure (Fig. [4^,) . 
These features may represent outflows from the central few 
parsecs of the galaxy. The bright Pa emission itself indi- 
cates the ionized inner portion of a dusty circumnuclear disk 
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around Sgr A*. This disk appears more extended, as traced 
by the 8 \im emission (Fig. [4^,). Detailed analysis/modeling 
of these structures is needed to further investigate the phys- 
ical state of the disk and the energetics of the GC cluster 
and/or the SMBH, which is currently in a radiatively ineffi- 
cient accretion state. 

Fig. [4fci presents a close-up of the Sickle HII region, re- 
vealing fingers of ionized gas resembling the "Pillars of Cre- 



ation" in M16 ( Cotera et al.|2006y . These fingers are appar- 
ently illuminated by ionizing photons from hot stars in the 
Quintuplet cluster ( Cotera et al.|2006 Simpson et al.|1997 



Cotera et al. in preparation). The Pa-emitting "Arches 
around the Arches cluster shows an amazing array of fine 
linear structures (Fig. [4J3) , which are preferentially oriented 
in the same directions as magnetic field vectors measured in 
infrared polarization ( |Chuss et al. 2003 Nishiyama et al. 
2009). Similar organized narrow linear Pa-emitting features 



are also abundant in some of the more compact HII regions 
seen in Fig. [2b. These structures may indicate a critical role 
that locally strong magnetic fields play in shaping the mor- 
phology and dynamics of the ionized gas (e.g 
et al |2009 \ 



Nishiyama 



Fig. [5 shows a selection of compact Pa-emitting neb- 
ulae. Some of these enigmatic nebulae may represent indi- 
vidual massive stars in formation (e.g., first row in Fig. (El), 
whereas others may be due to stellar ejecta and/or their 
interaction with dense ambient medium (second and third 
rows) . 



low up observations will at least double that number again. 
More importantly, over half of the newly identified stars are 
located outside the three known clusters. Therefore we have 
already more than doubled the number of known intercluster 
emission line stars. The newly discovered stars will enable 
us to address in greater detail the formation of massive stars 
outside clusters in our upcoming papers. In particular, the 
Pa star candidates outside the clusters are not distributed 
in a completely random fashion. In addition to the Pa emis- 
sion line, the F187N filter contains a few helium lines which 
are often strong emission lines in Wolf-Rayet stars. There- 
fore the most easily detected stars in the F187N filter are 
typically Wolf-Rayet stars, although future refinements to 
the data analysis will enable us to detect O and possible 
B supergiants as well, which have smaller Pa equivalent 
widths, in emission or absorption. Such refinements may be 
achieved by modeling the broad-band spectral energy distri- 
bution together with the narrow band measurements of the 
Pa emission and absorption. 



4 DISCUSSION 

While the data analysis, multi-wavelength comparison, 
follow-up observations, and modeling are still ongoing and 
will be presented in later papers, we discuss here the impli- 
cations of the survey for addressing several key issues related 
to star formation and its impact on the GC environment: 



3.2 Point- like sources 



4.1 How do massive stars form in the GC? 



The survey also enables us to perform an unprecedentedly 
uniform, high resolution, photometrically sensitive census 
of the stellar population. Our preliminary detection based 
on the F190N continuum image gives a total of about 0.6 
million point-like sources, which accounts for about 86% of 
the total observed F190N intensity in the survey field. This 
is the highest ratio of resolved point sources to continuum 
emission achieved to date in the very crowded GC field (e.g., 
Figs. [2^, and [3^,). The median of the unresolved intensity is 
~ 0.04 mjy arcsec -2 , averaged over the surveyed field. Our 
50% source detection limit is about 17th mag in the F190N 
band. A preliminary analysis shows that we should detect 
most of the massive main-sequence stars (earlier than B3, 
corresponding to ~ 8Mq) in the region as well as the bulk 
of red giants, including red clump stars. A detailed analysis 
of the spatial distribution of the detected sources and the 
still unresolved intensity will allow us to infer the overall 
intrinsic structure of the stellar light distribution in the GC 
and to provide quantitative measurements of the foreground 
extinction. 

Our preliminary analysis has resulted in the detection 
of about 200 Pa-emitting star candidates. Before the publi- 
cation of such candidates, however, we still need additional 
careful examination, especially for relatively weak sources, 
to minimize potential systematic contamination (e.g., due 
to local enhancements of diffuse Pa complexes). We easily 
detect all previously known emission line stars throughout 
the survey region (~ 70 total; mostly Wolf-Rayet stars). We 
have already verified an additional 30 emission line stars via 
spectroscopic follow-up observations ( |Mauerhan et al.|2009| 
Cotera et al. in preparation) . We anticipate that further fol- 



Two large questions remain unanswered about massive star 
formation in the GC: (1) Do massive stars form exclusively 
in the clusters, or in looser associations throughout the field? 
and (2) Is star formation continuous, or does it happen in 
discrete bursts? Until recently, massive (e.g., emission-line) 
stars have mostly be found within the 3 GC clusters. A num- 
ber of searches for young massive stars at infrared and X-ray 
wavelengths have been carried out (e.g., |Cotera et al.|1999 



Muno et al."p00 6a Mauc rhan, Muno fc M orris 2007| |MaueTl 
han et aL|2 009 ) . These studies, each with diverse search cri- 



teria, have yielded a rather incomplete picture of the young 
massive star population in this GC. 

The Pa emission is a sensitive tracer of massive stars 
and their environs. The emission traces warm ionized gas, 
produced and maintained primarily by Lyman continuum 
radiation from massive stars. Such gas can be associated 
with individual stars (e.g., massive stellar winds), their im- 
mediate vicinities (circumstellar material), and surround- 
ing ISM structures. The fact that the majority of our new 
sources are actually located outside the three known clusters 
has strong implications for understanding the dynamics of 
massive stars with ages of only a few x 10 yr. The exact time 
scales for clusters to be disrupted and the place where the 
stars end up provide information about the internal dynam- 
ics of the clusters as well as the large-scale tidal force of the 



GC ( 


Kim, Morris & Lee||1999| |Kim et al.||2000 |Kim, Figer 


& Morris 2004; Kim & Morris||2003| |Portegies Zwart et al. 


2002 


Portegies Zwart, McMillan & Gerhard 2003 


1. The spa- 



tial distribution of the massive stars detected in our survey 
can be used to better constrain these processes throughout 
the field. The distributions of the stars and their types as 
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well as follow-up work on both radial velocities and proper 
motions ( |Stolte et al.||2008| ) will allow us to constrain how 
quickly the massive stellar clusters are currently dissolving, 
how many star clusters have been disrupted, and whether 
or not the stars are formed in massive clusters which have 
been subsequently dispersed, or formed in relative isolation. 



4.2 What are the state and dynamics of the ISM? 

In the high pressure and high density environment of the 
GC, the diffuse Pa-emitting gas tends to be found at ion- 
ization fronts bordering molecular clouds. The gas density 
and pressure are expected to change drastically at such 



fronts (e.g., Scowen et al. 19981. The high resolution Pa 



map resolves various ionization fronts, which helps to deter- 
mine the local neutral gas density and the thermal pressure 
of diffuse warm ionized gas. The comparison between Figs. [I] 
and [2] demonstrates the complicated relationship among var- 
ious stellar and ISM components. Within the Sickle region 
(Fig. Bp), for instance, we see two distinctive ionization front 



morphologies. We see columns of photo-evaporating gas sim- 
ilar to the elephant trunks or pillars seen in regions such 
as M16 in one part, but also a remarkably smooth ioniza- 
tion front located on another side of the large HII region. 
Although very different, both are morphologically perpen- 
dicular to the brightest non-thermal radio filaments, which 
are uniquely found in the GC field and have been a long- 
standing mystery. They are known to trace enhanced mag- 
netic field and/or cosmic rays. But the underlying physical 
processes remain unclear. The observed relative morpholog- 
ical configuration then suggests that the large-scale, inter- 
cloud, vertical (poloidal) magnetic field is combing through 
the dense gas and gets enhanced in the process. The expan- 
sion of the magnetized inter-cloud material, probably accom- 
panied with cosmic ray acceleration, is likely driven by the 
strong mechanical energy output from the Quintuplet and 
Arches clusters, consistent with the barrel-shaped morphol- 
ogy of the nonthermal radio filaments (Fig. [I]). Therefore, 
the feedback from the clusters likely plays an important role 
in shaping the ISM and the eco-system of the GC in general. 
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Figure 5. Close-ups of Pa images of selected individual nebulae (equatorially projected). Detected sources have been subtracted, except 
for identified bright Pa emission star candidates, which help to show their potential relationship to the corresponding nebulae. The 
nebulae are labeled with rough Galactic coordinates, except for those with well-known names. The bar in each panel marks a 0.2 pc scale 
at the GC distance. 



4.3 Comparison with other observations 

A few Pa observations were made previously in the GC re- 
gion. The Sgr A West region, covering the inner ~ 4 pc of 
the Galaxy, was observed using both NIC3 and NIC2 (e.g., 
|Scoville et al.| [ 2003[ ). Several other programs concentrated 
only on the central parsec with NIC1. A NIC2 Pa observa- 



tion was also made for the core of the Arches cluster ( Figer et 
|al.|2002| ). The Pa emission from the Pistol Nebula (Fig.|4b) 
was covered with a 2 x 2 mosaic of NIC2 observations, which 



had a smaller field of view (19'.'2 on a side; Figer et al. 1999 1 



Our map shows considerably more structure in the diffuse 
Pa emission of the Nebula than the image shown in |Figer| 



et al. ( 1999 1, apparently due to the higher sensitivity of our 



NIC3 observations. In addition, our new data provide addi- 
tional epochs for variability analysis in these regions with 
the existing observations (e.g., Blum et al.|200l" i. 

Our study of the GC region complements other ef- 
forts to understand Galactic star forming regions, such as 



Orion (e.g., Feigelson et al.|[2~005 l, the Carina nebula (e.g., 



|Tapia et al.|2006| ), NGC 3603 (|Stolte et al.|2004[|2006l ), and 
Westerlund 1 (e.g., Muno et al.|2006b l. It is crucial to com- 



pare the time scales and efficiency with which stars form 
in the GC to those in other regions of the Galactic disk 
in order to gain a detailed understanding of the interplay 
between massive stars and the ISM. 

Ultimately, understanding the gas removal, via the for- 
mation of massive stars and their impact on the ISM or other 
processes (e.g., Type la supernova-driven Galactic bulge 
wind; |Tang et al^|2009| |Li, Wang fc Wakker|[2009), will 



help to determine why some host SMBHs accrete at rates 
near their Eddington limits, while others, including Sgr A*, 
are at rates orders of magnitude lower. 
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